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Abstract. We produce a hierarchiy of integrable equations by systematically 
adding terms to the Lax pair for the lattice modified KdV equation. The 
equations in the hierarchy are related to one aonother by recursion relations. 
These recursion relations are solved explicitly so that every equation in the 
hierarchy along with its Lax pair is known. 



1. Introduction 



The goal of this paper is to present an integrable hierarchy of equations associated 
with the lattice modified KdV equation (LMKdV). We do this by considering a 
Lax pair expansion which is polynomial in the spectral parameter. This kind of 
procedure is not new, in the same way hierarchies of integrable equations have been 
' produced [T]-[TS], including hierarchies of purely discrete equations starting with 

Where this work differs markedly from previously known results, is in that we 
can write formulae that express every equation in the hierarchy and its Lax pair 
•jH explicitly. 
^ In [17] two-component analogues of the LMKdV and the lattice sine-Gordon 

equation (LSG) were shown to represent the most general systems that can be 
associated with the Lax pairs for those equations. The two-component analogue of 
> LMKdV is 

. _ . . 

^ - + - = 3 + -, 1.1a 

\^ X y X y 

y X y X , , s 

^ + - = ^ + 1.1b 

y X y X 

where there are two discrete independent variables l,m G Z, the complex dependent 
variables are x = x{l, m), y = y(l, m) and shifts are represented by x = x{l -\- 1, m), 
^ X = x{l,m + 1) etc. This system was dubbed LMKdV2 because setting y = Xfi/x, 

where X{1) and fJ.{m) are arbitrary functions, retrieves the familiar LMKdV. A Lax 



pair for ( 1.1 ) is 



where v is the spectral parameter and the Lax matrices, L and M, appear in the 
linear system 

0^ Le, 9 = Me, 
which leads to the compatibility condition 

LM = ML. (1.3) 

In |18| an algorithm was developed to derive Lax pairs for the continuous mod- 
ified KdV hierarchy, based on the AKNS expansion technique [IH]- In this paper 
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we use a similar construction to form hierarchies of equations by adding terms to 



a generalisation of the L matrix in ( 1.2 1 as 



J _( ao + aiiy^ + . . . + aki'^'' boi^ + bn''^ + . . . + bki''^''+^ \ ,^ 
- [ coiy + c^,,^ + . . . + cW+' do + d,u'^ + . . . + dkiy^'' j' ^^-^^ 

where ai{l,m), <i < k are unknowns that wih be found using the compatibihty 
condition, other roman letters are similar. We keep the same generaUsed form of 
M for ah members of the hierarchy 

where the diagonal entries are set to unity using a gauge. 

While the form of the Lax pairs is nearly identical to that used in [TH], the 
methodology and results of this article are quite different. We will derive a set of 



recursion relations in the variable k from (1.4), which describe how the next Lax 
pair in the hierarchy can be found from the previous one. These recursion relations 
will be solved explicitly to give formulas for every Lax pair in the hierarchy as well 
as every equation. This will show that, as k is increased, the Lax pairs become 
increasingly complex but the equations retain the same simple form throughout. 
However, the number of arbitrary non-autonomous terms that can be inserted in 
to the reduced equations increases at each level of the hierarchy. 

Section [2] will see the compatibility condition put into a convenient form, then 
use this to derive recursion relations that relate different members of the hierarchy. 
In Section [3] we provide solutions to the recursion relations found in section |2] 
and describe how to use these to construct the Lax pairs including two examples. 
In section |4] we derive the associated nonlinear systems and show how these are 
reduced to scalar equations on stretched quad-graphs. The paper terminates with 
a discussion. 



2. Derivation of the recursion relations 

In this section we derive the recursion relations that relate different members of 
the hierarchy. In section [O] we put the compatibility condition into the required 
form and in section [2.21 we derive the recursion relations. 



2.1. Compatibility condition. Using a gauge to set ao ~ do ~ 1 and substituting 
(1.4) and (1.5) into the compatibility condition (1.3) gives one equation for each of 
the four entries of the 2x2 system: 

k 



(1,1) l + ^(a,+7&.-i)'^''+7&fc'^''+' 

i=l 
fe 

(2, 2) 1 + ^(d,; -I- /3c,_i)z.2* + Pcuv^''+^ 

i=l 

fe 

(1,2) Y.^pa^ + hV'+^ 
i=n 

k 

(2,1) ^(^d; + c,)i^2m 



l + Y,{a,+Pc,_^y+PckV^''+\ 

k 

i+Y,{d,+^b,_^y +^bkv^^+\ 

fe 

fe 

Y,{c,+-ia,)v^^+\ 



i=0 



i=0 



LMKDV HIERARCHIES 



3 



As the spectral parameter v is independent of the lattice variables I and to, we 
separate the compatibility conditions by powers of v. This yields 



Si + ^hi 


-1 




/3ci_i 


<i<k 


(2.1a) 


di + (ici 




= + 




0<i<k 


(2.1b) 


Pa, + 


k 


= Pd, 


+ h, 


<i<k 


(2.1c) 


idi + 




= la. 


+ Ci, 


0<i<k 


(2.1d) 



In addition to these we must consider the two extreme powers of v in the diagonal 
entries. The equations coming from diagonal entries of the compatibility condition 
at order i/^ are satisfied and at order !y2fc+2 permit the following parametrisation: 



h = y/x, 



Ck = x/y, 
1 = x/y, 



(2.2) 



where we have introduced the dependent variables x(l,m) and y(l,m), these will 
become the dependent variables in the nonlinear systems associated with the Lax 
pairs. 

Symmetry exists in the compatibility condition such that we may interchange 

{ [1 , 1] , [1 , 2] , [2, 1] , [2, 2] , a„ 6„ c„ , /?, 7} o { [2, 2] , [2, 1] , [1 , 2] , [1 , 1] , d„ c„ 6„ a„ 7, 

(2.3) 

where numbers in square brackets indicate the matrix entries of the compatibility 
condition, [1,1], [2,2], [1,2] and [2,1] are given by ( |2.1a[ ), ( |2.1b[ ), ( |2.1c[ ) and ( |2.1d[ ), 
respectively. Because everything in this paper is derived from the compatibility 
conditions, this symmetry naturally persists throughout. 

It is possible to remove one of the shifts in m fro m th e compatibility cond itions . 
This is done by starting at the i-th. condition from (2.1c), removing hi using (2.1a), 
Oi+i by using the i + 1-st in (2.1c), then successively removing the shifts in to by 
using the same equations. 



di 



-pd, + -h 



l_ 



1 



•i+i 



P 



'i+lj, 



etc. Because the last relevant compatibility condition, at order i/'^'^'^^, is 76/t — Pck, 
the procedure stops when i — k and we obtain 



P_ 
Pi 



E 



P 



Pc 



i+i 



(3d 



i+j 



Pi 



P 



An expression for hi is obtained in a similar way starting with (2.1al, we find 

k — i 



i"' f^.iPiV 



Pa 



i+j 



1 



H+j 



Expressions for Ci and di are obtained by symmetry (2.3). 

We write the compatibility conditions in matrix form by introducing 



(2.4) 
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which means that the compatibihty eonditions are equivalent to 

fc 



■0i 



1 ' 
1 



j^, + + V, , < * < fc, 
]=i+i 



(2.5a) 
(2.5b) 
(2.5c) 



where 
K = 



PIP 
ill 



l/(/37) 



0) 



l/(/37) 



= 

-1 PIP 
ill -1 

-7 7 



-^/(/?7) 
-7/(/?7) 1/7 



l/(/?7> 



J = 



^/7 
lIP 

-^/(/?7) 



G 



l/(/37) 



-7/(/?7) 1/7 

-1 ^/7 
7//3 -1 



2.2. Recursion relations. To solve the compatibility condition \2.b\ we perform 
the following sequence of calculations: 

• begin with the known quantity: tpQ — (1 1)"^, 

• use (2.5a)i=o to find (/Jq in terms of V'i and f i,i'>l, 

• use (2.5c)i=o to find V'l in terms of ipi, ipi and ipi, i>2, 

• repeat the las two steps at increasing values of i. 

Performing these calculations shows that solutions must be of the following form. 

fc 



(2.6a) 
(2.6b) 



for < i < fc, where Ai, Bi, P^, Q^/\ R!f \ S^^^ and T^^' become the unknowns for 
which w e^will derive re curre nce relations below. 

Take (2.66) and use (2.51 to remove ipi and (pi resulting in 



U), 



-E 



j=i+i 



Collect terms and note that 

j=i+lh=j + l j=i+2h=i+l 

to find 



(5['^^^X- J^f^^^-FV'- 
Bj - J - Ej=i+2 { 



)^,+i + (5['^^^i7 + rr^'j - jTr^' + Gf v.+i 



sy'H + f^^'^-jTy'-G. 



Ei=i+i ( 



sr'D 



(2.7) 
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Writing this in the same form as (2.6a)i+i and introducing 



leads to the following set of relations 



i-1 

E 

h=i+l 



s';'^H + Ty> - jry + G 



h=i+l 



(2.8a) 
(2.8b) 

(2.8c) 
(2.8d) 



The same argmiient can be followed starting with (2.6a) which results in the fol- 
lowing set of relations 



B, = A, ( A, - KA,j , 



= A 



0) 



A, 



(2.9a) 
,(2.9b) 



h=i+l 



(2.9c) 



where Ai = (^PiJ - KPi + . Equations J2.8i and I2.9I together form a set of 

nonlinear recursion relations that describe any Lax pair of the form (1.4) with an 



arbitrary number, k, of terms in each entry of the L matrix. If a solution to these 



recursion relations can be found, then tpi and (pi can be reconstructed from (2.6) 



which is a set of linear algebraic equations. In the next section solutions to the 
recursion relations are given. 



3. Solution to the recursion relations 

In this section we present exact solutions to the recursion relations derived above. 
These solutions enable us to construct every equation and Lax pair in the hierarchy 
explicitly (see section l4]). We use the solutions to construct the Lax pairs for 



the second and third members of the hierarchy in section 3.2 Some important 



alternative forms of the solutions are given in section 3.3 



3.1. The solutions. We introduce the following notation: 

a — a{l, m + h), 



In particular {l3j) = I3{1, m + h)j(l, m + h). 



Close inspection of (2.8) and (2.9) suggests the following theorem 
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Theorem 1. The set of recursion relations defined by {2.8) and {2.!^ are satisfied 
by the following formulae. 



A, 



B, = 



11^=1 P 7 



Let P, 



Pl.^ 

P2r 



we have 



j-l 2/1-2 



2/1-1 




3 = 1 

and P2S is the same except j3 and 7 are interchanged. 
We express the quantity Ti as 



2/1-1 • 



1 



2i+l 

1 - ^7/(/37) 



7 



with 



"~ 7 LV;?'2j+i 



J=i /3 7 - 

and r2,i is obtained from Fi ^ 6?/ interchanging fi and 7. 
Similarly = ^ ^ 



with 



Ai,.=/3n 



2^-2 2,-1 

J=i /3 7 



(3.1a) 
(3.1b) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



The matrices q[^\ R^^ , S'P^ andT'>^^ are all diagonal and the [2, 2] elements are 
obtained by interchanging /3 and 7 in the corresponding [1,1] elements. As such, 
we write only the [1,1] elements of these solutions below. The [1,1] element of s\'''^ 
is expressed as 



2i (h) 



(3.6) 



/i= 



(/3^)j 



where 



U. 



{h) 



dr(7) 



9=0 

/3n„=i 7 ^ 



(/37)' 



r2i 



n5=0,s#/i 
2i hi-l ''2(i-g)-l-l 



(/7) - m 



7 /3 . . . 7 (3 



EE- E 

/il=0/l2=0 'l2(i-s)=0 

(See notes on the definition of ^ after the end of this theorem.) 



(3.7) 
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The [1, 1] element of T^'''^ is expressed as 



2^ ^(h) 



(3.8) 



where 



and 



h=0 {P^)3 



3=0 



(3.9) 



/ll=0 /l2=0 



)— 'l2(i-g) + l— 



T/ie [1, 1] element of Q\''^ is expressed as 



2i-l (h) 

— 1i 



where 



^■'[1,1] h 
h=o 



g=0 

2i-l hi-l ''2(i-g)-2-l 'l2(i-g)-l 



(3.10) 



EE-- E E 'it.'-'r-' P 



hi=0h2=0 h 



2(i-g)-l-0 '«2(i~g) 

n<,=l 7 P 



The [1,1] element of is expressed 



l\.g=i3,g^h 

0) 



{h) - (A) 



h ' 



(3.11) 



where 



>) ^ 5:^(?^)(^)(/3^). 

5=0 

2i-l hi-l ''2(!-g)-3-l '»2(s-g)-2-l 

^^?^\(7) - E E--- E E 

/ll=0h2=0 /l2(i-g)-2=0 /l2(i-g)-l=0 



7 /3 . . . ^ 7 . 



Proof. The proof simply involves substituting the solutions into (2.8) and (2.9) to 
show that they are satisfied. However, it is rather technical and so it is left to 
Appendix 1^ □ 

Some notes on the definitions of S^^i^^ (7) : the upper index affects both the upper 
limit of the first sum and the number of factors (sums) , the argument of the function 
(in this case 7) refers to the first factor in the sum, factors alternate between j3 and 
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7 after that. The definition of il2i^\'^) is similar except there is one extra factor 
(and one extra sum). In general 



a hi—1 

EE- 

a hi — I 

= E E ■ 

/ll=0 ft,2=0 



1 p 



■ E E 

ha-b-l—^ ha-b—0 
ha-b-a — 'i- ha-b — 'i- 

■ E E 7/3- 

ha-b=0 /la-b+l=0 



p 7 



hi' 



where a — 6 > is assumed to be odd in the above formulae, but the only difference 
when a — & is even is that /3 and 7 must be swapped in the last two factors. 
Although the lower limits of all the sums are zero, we adopt the convention that 
12h=c f(^) ~ ^ when d < c, so the effective lower limit is greater than zero for 
many of the sums. For a = 6 we let ^i"' — 1, the case where a < b does not arise. 

3.2. Construction of the Lelx pairs. Here we explain how the solutions given 
in the previous section are used to construct the Lax pair for any member of the 
hierarchy. 



For fc = the Lax pair is already known, and is written immediately from (2.2) 



Starting with the formal solutions (2.6 1, it is obvious that the following linear 
system must hold when A; > 

Zk-^k - Cfc, (3.12) 

where 



Zk = 



I 


c.(l) 




o(2) 






1) 




1) 


q(k) ' 





/ 


Pi 






. Qt 


1) 




1) 










I 






. st 


1) 




1) 













I 


P2 ■ 


■ Qt 


1) 


H2 


1) 





and 







Pk-i 

I 




q{k) 



/ = 



1 
1 



^'fe = [ '/'O V'l 



1pk-l fk-l Ipk 



Ch — 



Bo 



t; 



(fe) 



(k) 



Vk 



Ak-i - R^k-iVk 



Bk-i - T[^\(pk 
Ak - Pk^Pk 

Both Uk and Ck are known from the formulae in the previous section so this system 



can be solved for ipi and (pi, 1 < i < k (i/jq and (fk are also known from (2.2)). The 
Lax pair for the fc-th member of the hierarchy is then obtained from ( 2.4 ) and ( 1.4 ). 
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As an example, consider fc = 1 where we have 



^1 



I S, 



(1) 



/ 

Using (|3JI), (13^, (|3J|) and (|3^ we find 



Ai - Piipi 



T 



— 

(1) 



-1/7 












1 

Pi 



1 
1 



Bo 



7/3 

/? 
7 



and 



Pi = 







(7/7+1)//? 

(/?/^ + l)/7 

Solving (3.12) and noting (2.2 1 gives us the terms in the L matrix from the Lax 
pair, which is 

1 + v^ixjx + xy/{xy) + x/x) viyjx + yjx + yjx) + v'^y/x 
iy{x/y + x/y + x/y) + v^x/y 1 + v'^{y/y + yx/{yx) + y/y) 
We also include the L matrix from the Lax pair when k = 2: 



L 



ao 



1, 



ai = 



0.2 



bo 



bi = 



X 


X 


X 


3 
X 


xy 


^3 

xy 
4! + 


:4 


X 


^3^ 
X 


^4 + 
X 


- + 
X 


3 ^ 

xy 




3 


X 

4 
X 


X 

1 1 

X 


xy 
xy 


4 

xy 

3 

xy 


3 

xa; 

+ — r 
XX 






y 




V 


3 


4 

2/ 






X 


X 


h 4- 

3 
X 


a; 


— 5 

X 






V 


V 


V 

h - -t- 


; 3 ^ 
xy xy 

.4 + " 
XX XX 


3 

2/2/ 


^4 
, 2/2/ 


3 
X 


X 


X 


4 i 
xy 


' 3 
xy 



3 

4 „ 
xy 



, 4 
xy 



4 
xy 



xy xy 



3 

xy 

XX 



2/2/ 



S 4 

_ 3^2/2/ 

_ ; + _ 3 > 

2^2/ xxy 



2/ 

X 



The other coefficients in the Lax pair, Ci and di, i — 0, 1, 2, are obtained from the 
corresponding bi and a^, respectively, by interchanging x and y. We then form L 
using (1.4). 

3.3. Alternative forms of the solutions. The main building blocks in the solu- 
tions are sums such as 



J=0 



(2i) 



2i hi-l 



J=0 



/ll=0 /l2=0 



• E 7/3 

'!-2(i-3) + l-0 



'»2(i-j) + l 

. /3 . 



In this section, we survey some properties of these sums. 
We make the claim 

Y.^'in)m^-Y. E ■•• E 7/^^ 

i=0 hi=0 11-2=0 /i2i+i=0 



h2i+i 

■ 7 , 



(3.13) 



where 5^ is the kroenecker S. Evidence to support this claim follows. It can be 
helpful to display all of the terms in such sums diagrammatically, an example of 
which with i = i, g — 5 is given in figure [ij where each non- increasing path from a 
node on the left of this figure, to one on the right represents a term in the sum. 



10 



MIKE HAY 




Figure 1 . Diagrammatic representation of sums such as those in 
(3.131. Any non-increasing path from a node on the left to one on 
the right represents a term in the sum. 



It is possible to separate the sums by removing terms containing the highest shift 
in m as follows 

hi=0 /i2i=0 hi=Q h2i-i=a 

2i^l h2,-i-l+s'^^'-^ 

til 



y h2 

hi=0 ft.2i=0 



E... E /3---7\ (3.14) 



or the lowest shift in m as follows 



/il=0 /l2i=0 /ll=l /!.2i-l=l 



2i /i2i-i-l+i5g^' 



E--- E /^■••7- (3.15) 

/ll = l /l2i = l 

Such decompositions are often useful and they appear again in appendix |Xj This 
process can be repeated and, when the number of factors reduces to two less than 
the difference between the limits of the first sum (the number of possible m-shifts), 
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a factor of {/Sj) can be taken outside the sum 

'i2i-l+'5o^' , L 2i h2i- 



2i 

E- 

hi = l 



-1+C 

E 

'i2i + l=l 



hi 
/3 



/3 =(/?7) 



2i 

E- 

hi = l 



-2-1+5/' 

E 

/l2.-l = l 



/3 



(3.16) 



To see that (3.13) holds, we proceed in this way until all of the repeated iterates 

a 

have been removed, producing a polynomial in {Pj) . Combining the coefficients of 
this polynomial to reproduce (3.131 is not a trivial exercise, but we omit it here 
because it is not essential to other parts of this article. 

A similar argument justifies the following expression for the sums involving ^: 



2i 

E- 

hi=0 



h2^- 



-1+S, 

E 



"1 h-, 
Pi 



(i 7 



(3.17) 



Equations (3.13) to (3.17) also lead to another alternative expression of the 
sums, this time as 2 x 2 matrix products, which can be useful for computer algebra 
calculations. Define 



hi=b 

where a — 6 is assumed to be even, when a 

swap from /3 to 7 or vice versa. If the argument of g), i = 0, 1, is changed 

from /3 to 7, then the first factor in the sums changes from /3 to 7 and the remaining 
factors alternate as usual. There are a — b factors (sums) in terms in uj^^l while 

uj^^l{f3;g) has a — b + 1 factors. When a = 2i and 6 = 0, (3.18) and (3.19) become 
the same quantities as those in (3.13) and (3.17) respectively. 
Using the same idea behind (3.14) and ( |3.16[ ), we find that 



E 



E 

ha-b+i=b 

h 

t-6-1— 1 + Og 



hi 



h2 



Pi 



ha-b+l 



ha-b=b 



An 2 

13 7 



(3.18) 



(3.19) 



b is odd the last factors of each sum 



V 1 7 

where a> g and b < g. Using the same idea behind (3.15) and (3.16) yields 




i(/3;5) 



i(/3;5) 



(3.20) 



(3.21) 



where a > g, b < g and a — 6 is even. By repeatedly applying (3.20) and (3.21) 
we can remove the sums in these quantities and replace them with 2x2 matrix 
products. Noting that 

( 4!](/3;5) 
V 4°ff(7;5) 

we find that 





(3.22) 



12 



MIKE HAY 



where products with a tilde, J^, are to be written in reverse order and 

if h — g is even, 













1 1 












i ) 



a h — g is odd. 



/3, if /i — 5 is even, 
7, if /i — 5 is odd. 



One can construct any of the relevant sums as matrix products by using (3.22) 



4. Associated nonlinear systems 



In this section we derive the nonlinear systems associated with the hierarchy 
of Lax pairs that are described by the formulae given in section [3] Because the 
compatibility conditions (2.1 ) from every order of the spectral parameter were used 



in the construction of the formulae of section [3] we can be sure that the entire 
system of compatibility conditions is consistent and we only need to check one of 
them (from each matrix row) to find the associated coupled nonlinear system. 



Noting (2.2), at i = k the [1,2] entry of the compatibility condition (2.1c) reads 

y - ' 



y 

X 



y 

X 



y 



(4.1) 



Using (2.6) and the formula (3.2), we find that 

flfe = 




k 

2k fi /-^ 

y " 3=1 



(4.2) 
(4.3) 



Substituting these values into (4.1) and into the corresponding compatibility con- 
dition from the [2, 1] entry (2.1d) gives us the coupled pair of nonlinear equations 
associated with this Lax pair. 



(4.4a) 
(4.4b) 



2k+l 
X 



y 



2k+l 

y 


X 




2k 

y 


X 


2k 

y 


2k+l 
X 


y 


X 


2k 
X 


= - H 

y 


~ 2k 
X 



In the case where k — Q, (4.4) is the LMKdV2 equation from [XT'. Indeed, (4.4) is 
simply LMKdV2 on a sheared lattice, as depicted in Figure [2] 

LMKdV2 is reduced to the familiar LMKdV by setting y = Xji/x where X{1) 
and l-i{m) are arbitrary functions. We can perform a similar operation in equations 
(4.4), setting y = a/x where a{l,m) must satisfy 



fc - _ fe+i 
a a = a a 
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/c = 



2fc+l 
X 



2k\ 

X I 




Arbitrary k 




Figure 2. Elementary quadralateral of LMKdV2 for fc = (left) 
and the points occupied on the sheared lattice by a member of the 
hierarchy with arbitrary k (right). 



so that each of the equations in (4.4) reduce to the same system. This implies that 
a must be of the form 

fe 

a[lm)=\(l)\{U{kl + mf"^ , (4.5) 

i=l 

where C,^ ^ 1 and /i, 1 < « < fc, as well as A(^), are arbitrary. Therefore, the 
reduced equations contain fc + 1 arbitrary functions. 

If we change variables such that I' = / and m' = kl + m then a'{l',m') = 
X'{l')fi'{m'), A' and fi' arbitrary, and exactly LMKdV is retrieved. 

5. Discussion 

We have constructed a hierarchy of integrable partial difference equations associ- 
ated with LMKdV. We have obtained explicit formulae that describe all equations 
in the hierarchy and all of their Lax pairs. We believe this is the first publication 
of such a result. 

It is remarkable that the equations retain the same simple form throughout 
the hierarchy, even though the Lax pairs become very complicated and contain a 
large number of terms as the order of the polynomial in the spectral parameter is 
increased. Indeed, the solutions to the recursion relations from section |3] contain 
far more terms than the Lax pairs, all of which file together perfectly to form the 
equations. 

At present the solutions to the recursion relations from section [3] are somewhat 
mysterious. The recursion relations are nonlinear and it will be interesting to 
understand how these solutions relate to the other known solution types. 

LMKdV is not a special case, there is evidence to suggest that it is possible to 
construct similar hierarchies for other lattice systems with 2x2 Lax pairs, which, for 
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example, includes all of the equations on the ABS list [ST] . Further work may also 
yield similar results in the widely studied continuous and semi-continuous domains. 
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Appendix A. Proof of Theorem [T] 



This appendix demonstrates why (3.6), (3.8), (3.10) and (3.11) are solutions to 



(2.8) and (2.9). Because it is the least complicated, we concentrate on the [1,2] 
entry of (2.8c), but the other parts are done in a similar way. 



Substituting (|3.6|) and (|3.8|) into (|2.8c|) returns 



( 



^i+l 



V 



9+1 , . 



9+1 , . 



Pi 

9 + 1 

13 7 



iP-ry 



0=0 l 3 + 1 



(9) 



9 + 1 • 

01- Pi 



Pl_ 
3 + 1 
Pi 



i-i 



-32 



(9) 



- 7t 



(9) 



(A.1) 

where we have set si'' = s^, while is obtained from si^ by interchanging /? and 
7, and follows the same pattern. Note that since we will conce ntra te on the [1,2] 
element of Q'f^i, the first column of the matrix on the RHS of (A.l) is irrelevant. 
To derive the above equation we have used 



j-i 

h=i+l 



1 



J-1 

h=i+i " 



2i 



-1 

1 



j-i 

E ffl^)j-'i 9+1 



PIP 
-1 

-Si 



(9) P^ (9) 



7 (9) 



P 



'32 



(9) 



2^ ^ 
" E g+1 

9=0 ^p^y 

and similar expressions where necessary. 



(9) P 



9+1 



5. (f) 

7''2i 



-^31 
-32 



(9) 

i 

(9) 



Expanding (A.l) and considering only the [1,2] element shows that 
1 „ ^ 1 „ 1 



= 



2i+l 

( /37 )■'■ 



2i 

E— 

9=1 {hY 



where 02i+i, 6g, < 5 < 2i, and Oq, which are given explicitly below, must all be 
identically zero because of the independence of j in this expression. We therefore 
have three equations which must be satisfied, these are 

9-1 /? 7 



e 



2i+l 



= /37s^ 



(2») 



2«+l ^f2il 
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+ -((/37)~/37)s2f^-(/37)ti. 
7 



1 * 

4n 



2/1-1 



2h 



P 7 



(9-1) 



- 7i; 



(9-l)x 



and 



2ft-l 



e, = o = In 



fl 2/1 (0) 

/3 7 S2- 



f, 1 fl' 2/1+1 



/37 



9+1 9+1 

9=0 [/3^-(/37)](/37)^ 



7(/37)* 

^ ^■Sli 



2/1-1 



n 

V h=l 



2h 



(A.3) 



(A.4) 



J 7 

2/i 



2/1+1 / \7 



/3 7 



1^(9) 7r(9) 

-S2i ~ -t2i 



First we show that (A. 2 1 is satisfied. Using (3.6) and (3.8), as weU as (3.13) and 



(3.7) we obtain 



/ll = l /l2i = l 

2i + l /i2i-l-l + 52i + l 



/ll''2 /l2,-l''2 

7 /? ■ 



7 /3-2^--- 2^ 1 P ■■■ 13 7 + 

/ll = l /i2i + l=0 

, 2i+l ''2! — l + l52^J:i , 

''J/12 %-'h2.^ 7 >r- 't- ''2 />2,'*^' + ^ 

P hx = l /l2. + l = l 



2i+l v-^ 

T E ■ 

/ll=l 



E 

/i2,=l 



However, the sums arising here can be re- written using (3.14), (3.15) and the ar- 
3 derive thos 

2i + l 'l2i-l + <52i + i 



guments used to derive those. Therefore, we may write the second sum from above 
as 



E-- E 7. 

'i2i + l = l 



7 



2i + l 'l2i-l-l + l52i + l , 

"1 /lo- 

7 



T E ■■■ E 



hx = l 



/!.2i = l 



and the fourth sum is similar. This shows that (A. 2) is satisfied. 



Turning to (A.3), we agian substitute (3.6) and (3.8), use (3.13) and (3.7), then 



cancel common terms to arrive at 

2i+l /i2i-i-l+5g^'"^ 



2i+l 



= - p lY. 
/ii=i 

„. , 1 2i+l /i2i-l+<5: 

- /3 E ■■■ E T- 



, 2j+l /i2i-l+i5, 

7.../3+7E--- E /3 

/l2i — 1 /li — 1 h2i-\-l—1 



hi /l2i + l 

/3 



/l2i + l 

7 



2i+l /i2i-i-l+i5g 



m E • 

hi—l /l2i+l — 1 /li — 1 

2^+1 9 ^ 

[(/37)-(/37)] E ■•• E /3...f. 

/il=0 /i2i=0 



E 

'l2i = l 



hi 



h2i 



/? . . . f 



It is easy to see that this cancels by using (3.14), (3.15) and (3.16). Thus, (A.3) 
is satisfied. We also note that equations (3.13) and (3.7) are not essential to show 
that (A.2) and (A.3) are satisfied. Although it is more complicated, we can use the 



solutions (3.6) and ( 3.8[ ) directly. Doing so leads to a cancelation of terms similar to 
that described in Figure [3j see (A.7), (A.8) and (A.9), below, and the surrounding 
arguments for an explanation of this figure. 
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Finally, consider (A.4|. Using the (3.6) and (3.8) we can rewrite (A.4| as 



2i 



-('^)fl(/37)&7)-E 

J = l 9=0 



9+1 



9+1 , 2i+l 



n;roWi[(/37)-(,97)] 

2i+l 



(A.5) 



In the following, we often need to differentiate between terms that contain mul- 
tiple factors shifted by the same amount in m and terms that do not. Let us refer 

9 

to any /3 and 7 shifted by the same amount in to, i.e. {P^), as 'repetitively shifted'. 
Let us refer to those terms that do not contain any /3 and 7 shifted by the same 
amount in to as 'cascading'. Both repetitively shifted and cascading terms can be 
seen in Figure [l] In that figure, terms containing repetitive shifts must all have at 
least two nodes along the horizontal line through g, while cascading terms have at 
most one node along this line. 



To show that (A.5) holds, we factor out all of the repetitively shifted terms and 



focus on the cascading terms. The strategy is to show that the cascading terms 
group together in certain ways and that the repetitively shifted terms multiplying 



these groups all cancel, proving that (A.5) is satisfied 



In (A.5), the quantities that contain only cascading terms are: 



Cr^(/3), Tr)r^(7), ^er^(7), Iv'^I-HPI 0<h<2^ 



(A.6) 

2i + l 2i+l 

(note that some terms within contain a factor of {l3j), which is taken 

into consideration below). 



Obviously we may use the definition (|3.9|) and write 
^^f^(7) 



^2(i-h)-l- 



2i hi-1 

EE- E 



-1 'l2 



E 7 ^ 



hi=0 h2=Q 



'2(i-h) 



=0 h. 



2i 



2i-l 



l2{i-h)' 



7E-- E 

/l2=0 'i2(i-h) + l=0 

2i-l '»2(i-)i)-l 



■ P 7 



P 7 



hi=0 



1=0 



{2h), X 

= (7) 



hi=2i 



P 7 , 



Q<h<2i. 



(A.7) 



hi<2i 



Notice that when h — 0, i}^^ has 2i + 1 sums, each with a decreasing upper limit 



Thus, ?72i' must have hi — 2i, so the second term in (A.7) is zero when h ~ 0. The 



other terms in (A.5) involving rj can also be decomposed in a similar manner. 
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We may decompose the terms involving ^ into parts that are related to rj as 
follows. 

2i '»2(i-fe)-l-l , h , . , 

1^ 1^ ■■■ 2^ P ^ ■■■ P 7, 

ftl=0/l2=0 'l2(i-h)=0 



2i-l 'i'2{i-h)-l- 



'l2(i- 



13 



1 



h2=0 



2i-l h2(i- 



?ll=0 'l2(i-fi)=0 



„. 2i-l 'i2(i-fi-l)-l 

'll=0 'l2(i-h-l) + l=0 



'l2(i-fe-l) ft, 

, /3 



2(i-h-l) + l 

7 



l2(i-fi)- 



(/?7) ''2=0 ''2(i-h) + l=0 

We may also decompose ^ in a different way as follows. 

1-1 . 



P 7 , 



E--- E E ^7 



h2(i-h)-l fi2 

■ P 7 



fti=0 'i2(i-;i)-i=0 /i2(i-h)=0 

7Z^--- /^7--- 

hi = l 'i2(i-fe)-l = l 



2i ''2(i- 



-E--- E ri 

hi = l 'i2(i-h)=l 



l2(i-h)-l h 



2(i-h) 



2i-l 

7E 



''2{i-h-l)~ 



E z'" 



'^■2(i-h-l) + l 



(A.8) 



-^7E- E >^ 

{P^) ''2=0 'i2(i-h) + l=0 

2i 



P 7 , 



-€'^^(7) 



hi<2i 



(A.9) 



/ii=2i 



Using (A.7|, (A.8) and ( |A.9[ ) it is easy to show that the cascading terms group 
together as shown in Figure |3l note that any repetitively shifted factors are omitted 
from that figure. We give directions about how to use this figure below. 

Before we do that, let us recall two so-called 'Euler Identities' which are crucial 
to what follows. For an arbitrary set of ? + 1 distinct points aj ^ C, < j < i, the 
first identity is: 

fl«.-E-f^^ \ „ (A-10) 

j=0 g=0 O-g llh=Q,h^gy'^h Agj 
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and the second identity is: 



g=0 U.h=0,h^gi^h - ag) 



= 1, 



(A.ll) 



for an integer p with Q < p < i. These two identities are either given explicitly in, 
or easily derived from arguments outlined in, section 2.7 of ^2] . 

Next we turn to the coefficients of each of the groups of cascading terms, all 
such coefhcients are made up of repetitively shifted terms. Consider the top level 
of Figure [3] which corresponds to /i = 0. This h also corresponds to the index from 
the s um in the third term, and the inner sum in the last term, of (A. 5 1. Due to 

2i+l 



(A. 



and (|A.9[) we find that £,fhl3) 



P "n^i (7) iP) have common terms 

in their decompositions, as depicted by a dotted line in Figure [3j The coefficients 
multiplying these terms in ( A.5[ ) sum to: 

, 2h 2/1+1 

n ^ , I m nLi (/37)(/?7) 

2i+l „. h 



m mum -Pi] 

, _ * 2h 2/1+1 

K/37) n mm 



2i 

E 



pi/m[i 



g+l 2i+l 

m/m] 



II m 



h=0 



h=l 



2i 

E 

g 



y.2^+1 

11/1=0, 

1 



-Am 



2g+l ■ 

m\ 



° m nto. 



\m " m\ 



(A.12) 



This is the same as (A. 10 1. 

The other parts of (A. 5), at the level of ft. = in Figure [s] are cancelled in a 
similar way. We turn our attention to the arbitrary level < ft < 2z of that figure. 
As an example, on the right hand side we observe the following set of terms with 
common factors: 

2i+l , 

^2i 



7 P ^ij). 



1V^'^'\P) 



iT\P)- 



hi<2i 



Adding together the relevant coefficients from (A. 5 1 leads to 



= 



Y{%Am~Pi] 



2 J 2j + l 



HPi)\[mm 

i=i 



J7+1 

Pi/m 



E- 



9+1 

m] 



[mt 



^'[mim-i] 



= 



2i+l 

E 



[mf 



{P%)] 



< ft < 2i. 



This is the same as ( A.ll ). One can check in a similar way that all of the parts of 
(A. 5) are either identically satisfied or equivalent to one of (A.lOl or (A. 111. We 
point out that the solutions match the truncation of the sums in (A. 5 1 because 

S2i ^ ^- 



Thus, the (1,2) entry of (2.8c) is satisfied by the solutions (3.6) and (3.8|. In a 
similar way, it can be shown that each of the recursion relations in (2.8) and (2.9) 
are also satisfied by (3.6), (3.8), (3.10) and (3.11). 
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h = 



2i+l 



-~>2i + l , . 

/3 €^(7) 



2i + l (r.. 



7??2?^(/3)ki<2i 



hi<2i 



2i+l 

7 /3 e^°H7) 



ft = 1 



2i+l , s 



2i+l , . 



2i+l , . 



P 




'lfl?hP)\hi<2t 



.-ivfi i.P)\hi=2t 



l€\P) 



ivif iP)\hi<2i 



■"■7??2i^(/3)k=2* 



\ "iflif iP)\hT<2iL 



2i+l 



7 P (7) 



(2i+l \P_j2i '{Vh,<2i 

P 4 (7)i7f?,,,, 2 



klFcS) 7€'^(/3)k.<2. ^ ^ 

2^+1, (sCn^ -vn^^'V/siL v-n^'^si 

P (7)k=2» /7V2i lPJ|/Ji=2i I ) 

' ■/ .{2h+2),j ''■■■^f|?^'\P)\h.=2^ 



^ ^^f^(7)h.=2^^> (7i.=2.] 

^2^ (7)fti = 2z.-- 77?^. ^ ^(^)|/,i<2i 

\7^2i (/3)k<2i 

7%i (/3)|fti<2i 



v J 



ft + 1 



hi<2i 

2i+l 



S'^'\P) 



2i+l 



/3 '}^f+''(7)k,=2,.../ -.(2ft+4) 



2i+l 



7Vk 



\p)k 



<2i 



hi<2i . 



\7^^f^^^(/3)|ft.^ 



7 /3 (7) 



h = i 



2i+l ,„.. 



2i+l 



'P f^^in) 

2i+l 

-p 



hi<2i 
hi=2i 



<2i 



'■■■7€'^(/5)k. 
iri^^\P)\h^=2i 



&\P) 




'S^{P) 


■ / 



2i+l 



7 P ^in) 



i4f\P) 



Figure 3. Relationship between the various terms that appear in 



(A.5). Terms that appear exphcitly in the equation are boxed. 
Sohd arrows indicate decomposition of a sum into two terms. 
Dashed Unes hnk groups of Uke terms. The variable h is the same 



as that from the inner sum of the last term in (A.5 1 
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